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ABSTRACT 

By cross-correlating AKARI all sky survey in 6 infrared (IR) bands (9, 18, 65, 90, 140, and 
C\| , 160/im) with the SDSS galaxies, we identified 2357 infrared galaxies with a spectroscopic 

■ redshift. This is not just one of the largest samples of local IR galaxies, but AKARI provides 

crucial FIR bands in accurately measuring galaxy SED across the peak of the dust emission 
at > lOOfim. By fitting modern IR SED models to the AKARI photometry, we measured the 
total infrared luminosity (Lm) of individual galaxies. 

Using this Lm, we constructed the luminosity functions of infrared galaxies at a median 
redshift of z=0.031. The LF agrees well with that at z=0.0082 (the RBGS), showing smooth 
and continuous evolution toward higher redshift LFs measured in the AKARI NEP deep field. 
By integrating the IR LF weighted by Lm, we measured the local cosmic IR luminosity 
density of fl IR = (3.8±^) x 10 8 L Mpc~ 3 . 

We separate galaxies into AGN (active galactic nuclei), star-forming, and composite by 
using the [Nil]/ Ha vs [OIII]/Hj3 line ratios. The fraction of AGN shows a continuous in- 
crease with increasing Lm. from 25% to 90% at 9< logLm < 12.5. The SFR# Q and L^oiin 
show good correlations with Lm for SFG (star-forming galaxies) and AGN, respectively. 
The self-absorption corrected Ha/H/3 ratio shows a weak increase with Lm with a sub- 
stantial scatter. When we separate IR LFs into contributions from AGN and star-forming 
galaxies (SFG), the AGN contribution becomes dominant at Lm > lO 11 ^©, coinciding 
the break of the both SFG and AGN IR LFs. At L IR < 10 n i Q , SFG dominates IR LFs. 
Only 1.1±0.1% of tt IR is produced by LIRG (L IR > 10 n L Q ), and only 0.03±0.01% is by 
ULIRG (Lm > 10 12 Lq) in the local Universe. Compared with high redshift results from 
the AKARI NEP deep survey, we observed a strong evolution of ^m G oc(l+z) 41±0 - 4 and 
Clf R N oc(l+z) 41±0 - 5 . Our results show all of our measured quantities (IR LFs, L* , ^m N , 
rtj^ G ) show smooth and steady increase from lower redshift (the RBGS) to higher redshift 
(the AKARI NEP deep survey). 

Key words: galaxies: evolution, galaxies:interactions, galaxies:starburst, galaxies:peculiar, 
galaxies:formation 
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Figure 2. Redshift distribution of the AKARI-SDSS galaxies. 
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Figure 1. Angular separation between the AKARI FIS and SDSS positions. 
We use galaxies matched within 20 arcsec for LFs. 



1 INTRODUCTION 

To understand the cosmic history of star formation and AGN (ac- 
tive galactic nuclei), it is vital to understand infrared (IR) emission; 
the more intense star formation, the more deeply it is embedded 
in the dust, hence, such star formation is not visible in UV but 
in the infrared. Similarly, AGN evolutionary scenarios predict that 
AGN are heavily obscured at their youngest, compton-thick stage 
dTreister. Urrv. & Viranill2009l) . The Spitzer and AKARI satellites 
revealed a great deal of infrared emission in the high-redshift Uni- 
verse, showing strong evolution in the infrared luminosity density 
(Le Floc'h et al. 2005; Perez-Gonzalez et al. 2005; Babbedge et al 
2006; ICaputi et alj 1200 71; Magnelli et al. 2009; Gruppio ni et al 



2010; Rodi ghiero et al.1 120101) . For example, at z=l, iGoto et al 
d2010d) estimated 90% of star formation activity is hidden by dust. 

However, to investigate evolution of IR LF, these high-redshift 
studies need a good comparis on sample at z=0. Even today, the 
often used is the IRAS LFs dRush. Malkan and Spinoglid 1 19931 : 
ISanders et af1l2003l; IGoto et alj|2010ch from 1980s, with only sev- 
eral hundred galaxies. In addition, for more than 25 years, bolo- 
metric infrared luminosities (Z/.r.R, 8-ioooum) of l ocal galaxies have 
been estimated using equation in IPeraultj dl987l) . which is a sim- 
ple polynomial, obtained assuming a simple blackbody and dust 
emissivity. Furthermore, the reddest filter of IRAS is 100/im, 
which does not span the peak of the dust emission for most 
galaxies, leaving a great deal of uncertainty. A number of studies 
found cold dust that c a nnot b e detected with the IRAS. For ex- 
ample, iDunne & Ealej feOOlh detected such cold dust with T ~ 
20K using Scuba 450,850 fim flux. ISvmeonidis eTai] j2009h de- 
tected cold galaxies with SED peaks at longer wavelengths using 
Spitzer/MIPS. These results cast further doubt on Lir estimated 
using only < 100/im photometry. More precise estimate of local 
Lir and thus the local IR luminosity function (LF) have been long 
awaited, to be better compared with high redshift work and to un- 
derstand where the end-point of the cosmic IR density evolution 
is. 

AKARI, the Japanese infrared satellite (Murakami et al. 

120071) . provides the first chance to rectify the situation since IRAS; 



AKARI performed an all-sky survey in two mid-infrared bands 
(centered on 9 and 18 fim) and in four far-infrared bands (65,90, 
140, and 160/im). Its 140 and 160/im sensitivities are especially 
important to cover across the peak of the dust emission, allowing 
us to accurately measure the Rayleigh- Jeans tail of the IR emission. 

Using deeper data and modern models, in this work, we aim 
to measure local Lir, and thereby the IR LF. By matching the 
AKARI IR sources with the SDSS galaxies, our sa mple contains 
~2357 IR galaxies. Co mpared to the previous work dSanders et al.1 
l2003l ; lGotoe"tal]|2010d) . the sample is several times larger, allow- 
ing us to accurately measure IR LF of the local Universe. The op- 
tical spectra of the SDSS also allows us to separate IR emission 
from AGN and star-forming galaxies. This work provides us an 
important local benchmark to base our evolution studies at high 
redshift both by the current AKARI, Spitzer, and Herschel satel- 
lites, and by next-gene ration IR sate l lites s uch as WISE, JWST, 
and SPICA. Following ISanders et alj d2003l) . we adopt a cosmol- 
ogy with (h, fi m ,fi A ) = (0.75, 0.3, 0.7). 



2 DATA AND ANALYSIS 
2.1 AKARI-SDSS sample 

AKARI iMurak ami et alj [2007) performed an all-sky survey in 
two mid-infrared bands (centered on 9 and 18 /im) and in 
four far-infrared bands (65,90,140, and 160/im). In this work, 
we use the version 1 of the AKARI/IRC point source cata- 
log and the AKARI/FIS bright source catalog, which is se- 
lected in 90/im. The 5 a sensitivities in the AKARI IR filters 
(SW, L18W, N60, WS, WL and 7V160) are 0.05, 0.09, 2.4 , 
0.55, 1.4, and 6.3 Jy dlshihara et aitold ; lYamamura et al.ll2009l) . 
In addition to the much improved sensitivity and spatial resolution 
over its precursor (the IRAS all-sky survey), the presence of 140 
and 160/im bands is crucial to measure the peak of the dust emis- 
sion in the FIR wavelength, and thus the Lir of galaxies. 

To measure Lir, we need spectroscopic redshift of individual 
galaxies. We have cross-correlated the AK ARI FIS bright sourc e 
catalog with the SDSS DR7 galaxy catalog dAbazaiian et alj 2009). 
which is the largerst redshift survey in the local Universe to date. 

We used a matching radius of 20 arcsec, as shown in Fig[TJ 
This radius is determined based mainly on the positional accuracy 
of the AKARI FIR sources since for the SDSS, positional accu - 
racy is less than 0.1" to the survey limit (r ~22 |Pier et alll2003l) . 
AKARI's PSF size in the 90/im filter is 39±1". The astrometric 
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accuracy of the AKARI FIS bright source catalog is 3.8" in cross - 
scan direction and 4.8" in in-scan direction dYamamura et all 2009). 
The 20" radius also accounts for the physiscal shift that could be 
present between the peaks of the FIR and the optical emission. 
At the median redshift of z=0.031, 20" corresponds to 12 kpc. 
The so urce density of the SPSS main spectroscopic galaxies is 92 
deg -2 dStrauss et al.ll2002h . Therefore, within the 20" around each 
AKARI sources, 0.009 of background galaxies are expected to be 
present, i.e., less than 1% of our sources are expected to be a chance 
coincidence, which should no t affect our estimate of the LFs. See 
dSufherland & Saunderslll992l) for a maximum likelihood approach 
to this. 

We do not use galaxies at z < 0.01 since (i) at such a low 
velocity, uncertainty from peculiar velocity is too large to estimate 
distance from spectroscopic redshift, and (ii) the AKARI's photom- 
etry of largely extended objects is not finalized yet. We also remove 
galaxies at z ^ 0.3 to avoid too much evolutionary effect within 
the redshift range. We also use SDSS galaxies with extinction cor- 
rected r pe tro ^ 17.7, where the SDSS main galaxy spectroscopic 
target selection can be considered complete. To obtain reliable LFs, 
we only use AKARI sources with Sgo M m > 0.7Jy, where com- 
pleteness is greater than 80%. Among these, 99% of sources have 
high quality flag of f qual=3 in 90/im. 67% of AKARI sources 
with Sgo/jm > 0.7 Jy in the area covered by the SDSS had an op- 
tical counterpart. The rest of the sources are either stars or galaxies 
fainter than r pe tro = 17.7, which are property accounted for the 
LF computation in Section [3~T1 

Previous sample of local IR LFs was based on IRAS Sao^m 
flux, and thus, it has been a concern that colder galaxies, whose 
FIR emission peak at larger wavelengths, might have been missed 
in the sample. Our sample is selected using Sgo^m flux, and thus, 
more robust to such a selection bias. 

We found cross-matches of 2357 galaxies, with which we will 
measure IR LFs. Among these, 97, 98, and 87% of sources have 
measured flux in 65, 140 and 160^im. For these source the AKARI 
IRC catalog is also cross-matched. Only 0.8% of FIS-SDSS sources 
do not have an IRC counterpart within 20 arcsec. Fig|2] shows a 
redshift distribution of the sample. The number o f galaxies used are 
similar to the IRAS based work bv lGotol l l2005l) . despite the twice 
increase in the optical (SDSS) sample. This is due to more stringent 
cross-matching criterion we used. AKARI's improved astrometry 
allowed us to remove ambiguous cross-ids. Nevertheless, thanks 
to the large sky coverage of the SDSS and AKARI surveys, this 
is one of the largest number of galaxies used to construct IR LFs. 
For example, previous local IR LFs a re based on ~600 galaxies 
dSanders et al.l2003l:lGoto et alfcOlOch . 



2.2 Estimating total IR luminosity 

For these galaxies, we estimated total IR luminosities (Lib) by 
fitting the AKARI photometry with SED templates. We used the 
LePhare cod43 to fit the infrared part (>7 ' fim) of the SED. 
We fit our AKARI FIR photometry with the SED templates 
from ICharv & Elbazl (CHEL hereafter I200H) . which showed most 
promising results among SED models tested bv lGoto et aLlfeOlQch . 
It is a concern that the CHEL models do not include SEDs for 
AGN. However, we show later in Section 2.3 that this is not a 
major problem using observed data. We did not include IRAS 
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Figure 3. Lib measured by the AKARI as a function of spectroscopic 
redshift. 

photometry in the SED fit partly because the better spatial res- 
olution of AKARI often resolves ba ckground cirrus and nearby 
companion better jjeong et"al] l2007t) . and partly because not all 
of the AKARI sources are detected with the IRAS due to the 
difference in flux limits. Although the shape of these SEDs are 
luminosity-dependent, the large baseline from AKARI observa- 
tions (S9W, L18W, N60, WS, WL and AT 160) allows us to adopt 
a free scaling to obtain the best SED fit, which is then rescaled to 
derive Lib,. To be precise, we only used 65/im flux to free scale 
the SEDs since mid-IR bands can sometimes be affected by the stel- 
lar emission. For those sources with low quality flag (f qual=l), 
we adopted a minimum error of 25%. In this work, Lib is mea- 
sured in the wavelength range of 8-1000/im. Fig.[3]shows Lib as 
a function of spectroscopic redshift. 



2.3 AGN/SFG separation 

One has to be careful that two different physical mechanisms con- 
tribute to total infrared emission of galaxies: one by star-formation 
and the other by AGN. It is therefore fundamental to separate IR 
contribution from these two different physics, to understand the 
evolution of the star-formation and AGN activity. However, this 
has been notoriously difficult and has been a subject of active de- 
bate over the last decades. Many AGN/SF separation methods have 
been proposed such as X-ray, radio luminosity, optical line ra- 
tios, PAH strengths, submm properties and so on dZhu et alj2008t 
ISaiina et alj|2008l;lLaMassa et alfeOld) . Often, not all of the indi- 
cators do not agree one another, partly due to complicated nature 
of IR sources (multiple core, composite ...etc). 

In this work, since we have optical spectra of individual galax- 
ies, we use [Nil]/ Ha vs [OI 1 1]/ H j3 line ratios to classify galax- 
ies into AGN or SFG (star-forming galaxies). In Fig|4] we plot 
[Nil]/ Ha against [OIII]/H f3 line ratios. The black contours 
show the distribution of all emission line galaxies in the SDSS, 
regardless of the AKARI detection. The dotted and dashed lines 
are th e A GN/SFG separation crite ria presented by iKewlev et al.l 
d200ll) and lKauffmann et alj fc003b . respectively. We regard galax- 
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Figure 4. Emission line ratios used to select AGNs from our sample. The 
contour shows distribution of all galaxies in the SDSS with r < 17.77 (re- 
gardless of IR detectio n). The dotted line is the criterion between starbursts 
an d AGNs described inlKewlev et aljfiOQll) . The dashed line is the criterion 
bv lKauffmann et al] HopK Galaxies with line ratios higher than the dotted 
line are regarded as AGNs. Galaxies below the dashed line are regarded as 
star-forming. Galaxies between the dashed and dotted lines are regarded as 
composites. The blue and red dots are for ULIRGs, LIRGs, respectively. 
The green squares are Seyfert 1 galaxies identified by eye-balling optical 
spectra. 



ies above the dotted line as AGN, below the dashed line as SFG, 
and those in between as composite galaxies. 

Although these line diagnosis criteria work well for narrow- 
line AGN, they could be problematic for broad-line AGN. Al- 
though, in principle, one could separate narrow line components 
from broad-line by multiple-line fitting, it is often difficult with 
weak lines on noisy spectra. To see how well broad line AGN are 
classified in Fig|4] one of us (M.M.) has eyeballed all emission line 
IR galaxies with FWHM(Ha) greater than FWHM([OIII]) by more 
than 150 km/s. We identified 22 broad-line AGN, which are shown 
in green squares in Fig|4] As expected, the broad-line AGN are 
scattered everywhere in the Figure and escapes from the above cri- 
teria. Therefore, we manually overruled the classification of these 
broad-line AGN, and treated them as AGN in the rest of the paper, 
regardless of the line ratios in Fig|4] In this process, we also came 
across 20 narrow line AGN (although these are not a complete sam- 
ple). They were successfully classified as AGN in Fig|4] 

In FigO ULIRGs and LIRGs are marked with red dia- 
monds and blue dots, respectively. It is interesting that major- 
ity of (U)LIRGs are aligned along the AGN branch of the dia- 
gram, implying the AGN fraction is high among (U)LIRGs. This 
is more clearly seen in Figj5] where we plot fractions of AGN 
as a function of Lir. The red solid line shows fractions of AGN 
only, and the black dotted line includes composite galaxies. Note 
that the fractions do not include completely obscured sources in 
optical for both SFG and AGN. The fractions of AGN (includ- 
ing composites) increase from 15% (25%) at logL_r_R=9.0 to 55% 
(90%) at logL/ji= 12.3. This r e sults agree with previous AGN 
fraction estimates dGotol 120051 : lYuan. Kewlev, & Sandersl l20ld : 
lKartaltepeetalJl2010h . Improvement in this work is that due to 



Figure 5. Fractions of A GN and compo s ite ga laxies as a function of Ljr. 
AGN are classified using iKewlev et all feOOll) among galaxies with all 4 
lines measured. Compo site galaxies include those classified as AGN using 
iKauffmann et alj 120031) . 
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Figure 6. SEDs of AGN (red triangles) and SFG (blue squares). The red 
and blue solid lines connect median values for AGN and SFG, respectively. 
Note that only galaxies detected in all the 6 AKARI bands are shown here. 
The flux are normalized at 90/^m. 



much larger statistics, we were able to show fractions of AGN in 
much finer luminosity bins, more accurately quantifying the in- 
crease. Especially, a sudden increase of /agn at logLm >11.3 
is notable due to the increased statistics in this work. 

Having AGN/SFG classified, in Fig[6j we show restframe 
SEDs of AGN (red triangles) and SFG (blue squares), separately. 
Note that only galaxies detected in all the 6 AKARI bands are 
shown here. The flux are normalized at 90^im. The solid lines con- 
nect median points for each sample. In Fig|6] there is no significant 
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Figure 7. The SFR computed from the Ha emission line is plotted against 
the total infrared luminosity (Lir). The Ha lines are corrected for stellar 
absorption and extinction using the Balmer decrement. AGNs are excluded 
from this figure using Fig(4] 
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Figure 8. The self-absorption corrected Ha/H/3 ratios are shown against 
Ljr for SFG. The solid line connects median values. 



difference between AGN and SFG SEDs. This has two implica- 
tions: although the CHEL SED models do not include SEDs for 
AGNs, the SED fit can be performed using the galaxy SED mod- 
els. At the same time, we need to be careful that although we se- 
lected AGN based on the optical line ratios, infrared SEDs of these 
optically-selected AGN can be dominated by star-formation activ- 
ity in the AGN host galaxy. 

In Fig[7] we show star formation rate (SFR) co mputed from 
the self-absorption and extinction corrected Ha flux dGotoll2005h 
against Lir for SFGs. The SFR is also corrected for the fiber loss 
of the SDSS spectrograph. Th ere is a correlatio n, but compared 
with the well-known relation bv lKennicuttl |l998), there is a signif- 
icant o ffset. This might warn us a sim ple applicat i on of iKennicuttl 
( 1998)'s law to IR luminous galaxies. Kennicutt 1 1998)'s sample 
was dominated by more regular (less IR luminous) spiral galaxies, 
whose Lir might stem from cool cirrus clouds warmed by older 
stars. This might not be the case for vigorously star-forming LIRG 
and ULIRGs with hot dust. 

In Fig(8] we plot the self-absorption corrected Hal Hf3 ratio 
as a function of Lir for SFG. Although the median in the solid line 
shows a slight increase with increasing Hal Hf3 ratio, but with a 
large scatter, suggesting that not only Lir decides the dust extinc- 
tion in galaxies. 

In Fig(9]we plot Lir vs Lyoiu}- Only objects classified as 
AGN are shown here. Both of these luminosities are considered 
to be good indicators of the AGN power, and there is a recog- 
nizable correlation, however, with a significant scatter. The exact 
reason of the scatter is unknown but possible causes include con- 
tamination from HII regions, cirrus emission, shocks, and emission 
from very little narrow-line region directly powered by an accreting 
black hole. 

In Fig[l0] w e show restframe u — r vs M r color-magnitude 
diagram. We used B la n ton et al.l {2003, v4_l) for k-correction. The 
contours show the distribution of all SDSS galaxies with r < 17.7 
regardless of IR detection. ULIRGs and LIRGs do not follow nei- 
ther the red-sequence at u — r ~ 2.6 nor blue cloud at u — r < 2.0, 
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Figure 9. Luminosity of the [OIII](5008A) emission line is plotted against 
total infrared luminosity for AGNs. Only objects classified as AGN are 
shown here. 

but widely distributed in the u — r color. As expected, ULIRGs are 
more luminous in M r than LIRGs. 

The red squares, green diamonds, and blue dots are for AGN, 
composite, and SF IR galaxies classified using Fig|4] These IR 
galaxies also show wide distribution in Fig[l0] with AGN slightly 
brighter than the composite, reflecting the fact that the AGN frac- 
tions are higher for more luminous objects (Fig|5]. The SFG are 
slightly more extended to bluer clouds than AGN and the compos- 
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Figure 10. Restframe u-rvs Mr color-magnitude diagram. The contours 
show the distribution of all SDSS galaxies with r < 17.7. The orange 
stars and crosses denote ULIRGs and LIRGs in our sample. The red square, 
green diamond, and blue dots are for AGN, Composite, and SF IR galaxies 
classified using Fig|4] 



ite. These distributions are consistent with AGN and IR galaxies 
populate the so-calle d "green valley" between re d-sequence and 
blue clouds JSmolcid2009l : ISchawinski et al]20ld) . Perhaps phys- 
ical reasons deciding u — r color are different for AGNs and SFG. 
For example, dust-reddening might play a significant role for SFG, 
especially infrared luminous ones. For AGN, lack of gas in red- 
sequence galaxies, and dominance of HII regions in blue-cloud 
galaxies may push AGN into the green valley. 

We caution readers that our AGN/SFG classification is by no 
means complete, but merely one attempt. In reality most of galaxies 
have both AGN and star-formation, and thus contribution to Lir 
and Lyoin] from both components. When one is dominant, it is 
difficult to assess a contribution from the other, which inevitably 
becomes a contamination to the former. Recently, it has been re- 
ported that to fit FIR SED of AGN dominant gal axies, one always 
need a starburst componen t , even for quasar s (IShao et al.l 1201 (j ; 
lHatziminaoglou et aljhoifl IPozzi et al]T2010l) . In this sense, our 
Lir of AGN include IR emission from an AGN host galaxy. The 
latter might be dominant in FIR wavelength. Those galaxies with 
line ratios consistent to be LINERs are treated as AGN in this pa- 
per, but their emission may have nothing to do with an accreting 
black hole, instead may be from shocks or something else. There- 
fore, it is possible that these classified as AGN in this paper still 
have significant infrared emission from host galaxies, instead of 
central black holes. To make the matter worse, it is additional un- 
known how many completely obscured galaxies we miss in our op- 
tical classification for both AGN and SFG. While this work had 
only optical line ratios for classification, the final conclusion must 
be based on multi-wavelength data. 



3 INFRARED LUMINOSITY FUNCTIONS 
3.1 The W m ax method 

With accurately measured Lir, we are ready to construct IR LFs. 
Since our sample is flux-limited at r = 17.7 and 5go M m = 
0.7 Jy, we need to correct for a volume effe ct to compute LFs. 
We used the 1/V max method dSchmidtH "l968). An advantage of 
the 1/Vinax method is that it allows us to compute a LF di- 
rectly from data, with no parameter dependence or a model as- 
sumption. A drawback is that it assumes a homogeneous galaxy 
distribution and thus is vulnerable to local over-/under-densities 
dTakeuchi. Yoshikawa. & Ishi i 2000). 

A comoving volume associated to any source of a given lu- 
minosity is defined as Vmax = Vz max — V Zmin , where z m in is the 
lower limit of the redshift and z max is the maximum redshift at 
which the object could be seen given the flux limit of the survey. 
In this work, we set 2 m i n =0.01 since at a very small redshift, an 
error in redshift measurement is dominated by a peculiar motion, 
and thus, Lir also has a large error. 

F or the infrared detectio n limit, we used the same SED tem- 
plates dCharv & Elbazll200lh as we used to compute Lir for k- 
correction to obtain the maximum observable reds hift from the 
Sgo U m flux limit. For optical detection limit, we used lBlanton et al] 
d2003l v4_l) for k-correction. We adoped a redshift where galaxies 
become fainter than eithter r- of Sgo^m-flux limit as z max . 

For each luminosity bin then, the LF is derived as 

1 1 

^ = al1^v — ~ Wl ' (1) 

L-^1j vmax, 1 

i 

,where Vmax is a comoving volume over which the ith galaxy 
could be observed, and AL is the size of the luminosity bin, and 
Wi the completeness correction factor of the ith galaxy. Note that 
the Vmax method considers both the optical and IR flux limits. 
Therefore, even an IR source with r > 17.7, a detectable vol- 
ume of such sources is properly considered by the Vmax method 
as long as IR sources with the same optical- to-IR ratio can be de- 
tect ed at a lower redshift. See Appendix of ISerieant et al.l d200lh 
and lAvni & Bahcalil (1980) for more details of a multivariate Vmax 
method. A very obscured IR galaxy (with r > 17.7 even at the low- 
est redshift of z=0.01), however, would have a detectable volume of 
zero, and thus, would not be included in the LF computation. The 
number density of such obscured IR galaxies is a matter of inter- 
est itself, but requires deeper optical data. Although we cannot test 
deeper r-band criterion since the SDSS spectroscopic survey stops 
the current limit, we can test a brighter r-band criterion. We com- 
puted the LF using a brighter r-band criterion by one magnitude, at 
r < 16.7. The resulting LF was noisier due to the reduced number 
of galaxies used, but essentially the same as the LF presented in the 
following sections. The test suggests that the Vmax method works 
well. We used completeness correction measured in the FIS bright 
source catalog release note for Sgo/im flux. This correction is 25% 
at maximum, since we only use the sample where the complete- 
ness is greater than 80%. Completeness correction in terms of sky 
coverage of both surveys is taken into account. 



3.2 Monte Carlo simulation 

Uncertainties in the LF values stem from various factors such as the 
finite numbers of sources in each luminosity bin, the k-correction 
uncertainties, and the flux errors. To compute these errors we per- 
formed Monte Carlo simulations by creating 150 simulated cata- 
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Table 1. Best double power-law fit parameters for the AKARI LFs. 



AKARI NEPD 



0.8<z<1.2 
0.5<z<0.8 
0.2<z<0.5 



AKARI SDSS (0.01 SzS0.30) 
AKARI RBGS (z = 0.0082) 



10.0 10.5 1 1.0 1 1.5 12.0 12.5 13.0 

log MU)] 



Sample 


L *IR ( L ©) 1 


)*(Mpc~ 3 dex _1 ) a (faint-end) /? (bright-end) 


Total 


5.7±0.2 x 10 1() 


0.00062±0.00002 1.99±0.09 3.54±0.09 


SFG 
AGN 


5.0±0.2 x 10 10 
4.3±0.2 x 10 10 


0.00043±0.00003 1.8±0.1 3.5±0.3 
0.0004±0.0001 1.4±0.1 3.1±0.2 


Table 2. Local IR luminosity densities by all galaxies, ULIRGs and LIRGs. 


Sample 


Qjotal (i Q Mp C - 3 ) 


Qf^iL^Mpc- 3 ) nY^ mG (L Q Mpc- 3 ) 


Total 
SFG 
AGN 


x 10 8 
5.9±{ t x 10 7 
l.9t° i x 10 7 


4.2t°; 7 x 10 6 L2±l\l x 10 5 
2.ll°; 7 x 10 6 UIqJ x 10 5 
1.9±°;| x 10 6 2.5±J;| x 10 5 



Figure 11. Infrared luminosity function of AKARI-SDSS galaxies. The 
L IR is measured using the AKARI 9,18,65,90,140 and 160/xm fluxes 
through an SED fit. Errors are computed using 150 Monte Carlo simula- 
tions, added by Poisson error. The dotted lines show the best-fit double- 
power law. The green dotted lines show IR LF at z=0.0082 by iGoto et al J 
<2010d) . The dashed-d otted lines are hi gher redshift results from the 
AKARI NEP deep field IGoto et all20103) . 
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Figure 12. An infrared luminosity histogram of galaxies used to compute 

FigE] 



logs, where each catalog contains the same number of sources, but 
we assign each source new fluxes following a Gaussian distribution 
centered at fluxes with a width of a measured error. Then we mea- 
sured errors of each bin of the LF based from the variation in the 
150 simulations. These estimated errors are added in quadrature to 
the Poisson errors in each LF bin. 



3.3 IR luminosity function 

In Fig[TT] we show infrared LF of the AKARI-SDSS galaxies, fol- 
lowed by Figl 121 which shows number of galaxies used to compute 
the LF. The median redshift of our sample galaxies is z=0.031. 
Overplotted in green dotted lines is IR LF of RBGS at z=0.0082 
but with re-measured L m based on the AKARI 6-band photom- 
etry jGoto et aL 201 Oct). Our LF agrees very we ll with that from 
the RBGS ISanders et al.l2003llGoto et alj2010d) . Although it was 
a concern that the RBGS was S^oum selected, the agreement sug- 
gests the selection does not affect the total IR LF, perhaps because 



is close to the peak of the dust emission and captures most IR 
galaxies. Dotted lines show higher redshift results from the AKARI 
NEP Deep field jGoto et alj2 10a). Our results (open circles) agree 
with that from z=0.0082, and show smooth evolution toward higher 
redshift. The faint-end slope of our sample agrees well with imme- 
diately higher/lo wer redshift LFs. 

Following ISanders et al] d2003l) . we fit an analytical func- 
tion to the LFs. In the literature, IR LF s were fit better b y 
a double-power law ( B abb edge et al l |2006|; IGoto et al] l2010allcj) 
or a double-exponen tial ( Pozzi etalT 20041 " iTakeuchi et al] 200( ; 
Lc Floc'h et al.ll2005l) than a Schechter function, which steeply de- 
clines at the high luminosity and underestimates the number of 
bright gala xies. In this work, we fit the IR LFs using a double- 
power law ( B abb edge et~alll2006h as shown below. 



$(L)dL/L* = $* 



${L)dL/L* = $* 



dL/L*, (L < L* 



1-/3 



dL/L* , (L > L*) 



(2) 



(3) 



Free parameters are L* (characteristic luminosity, Lq), <f>* (nor- 
malization, Mpc -3 ), a, and /3 (faint and bright end slopes), respec- 
tively. The best-fit values are summarized in TableQ] In FigQj] the 
dashed line shows the best-fit double-power law. The local LF has 
a break at i*=5.7±0.2 x 10 10 Lq. Understanding how this break 
(L*) evolves as function of cosmic time, and what causes the break 
is fundamental to galaxy evolution studies. This work provides an 
important benchmark in the local Universe. 



3.4 Bolometric IR luminosity density based on the IR LF 

One of the primary purposes in computing IR LFs is to estimate 
the IR luminosity density, which in turn is a good estimator of the 
dust-hidden cosmic star formation density dKennicuta l 1998). once 
the AGN contribution is removed. The bolometric IR luminosity of 
a galaxy is produced by thermal emission of its interstellar matter. 
In SF galaxies, the UV radiation produced by young stars heats 
the interstellar dust, and the reprocessed light is emitted in the IR. 
For this reason, in star-forming galaxies (SFG), the bolometric IR 
luminosity is a good estimator of the current SFR (star formation 
rate) of the galaxy. 

Once we measured the LF, we can estimate the total infrared 
luminosity density by integrating the LF, weighted by the lumi- 
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Figure 13. IR LF is separated for star-forming galaxies (blue triangles) 
and AGN (red squares), using classification in Fig|4] Total IR LF in Fig l 1 1 1 
is shown with the black dotted lines. 



nosity. We used the best-fit double-power law to integrate out- 
side the luminosity range in which we have data, to obtain esti- 
mates of the total infrared luminosity density, Qir. Note that out- 
side of the luminosity range we have data (Lir > 10 1 ' a L or 
Lir < 1O 9 ' 8 L ), the LFs are merely an extrapolation and thus 
uncertain. 

The resulting total luminosity density is Qir= (3.8jlf '%) x 10 8 
L Mpc~ 3 . Errors are estimated by varying the fit within la of 
uncertainty in LFs. Out of Qir, 1.1±0.1% is produced by LIRG 
(Lir > 10 n L Q ), and only 0.03±0.01% is by ULIRG (L IR > 
10 12 Lq). Although these fractions are larger than z=0.0081, still a 
small fraction of Qir is produced by luminous infrared galaxies at 
z=0.031, in contrast to high-redshift Universe. We will discuss the 
evolution of Qir in Section |4~Tl 



3.5 AGN/SFG IR LF 

Once we separated individual galaxies into AGN/SFG, we can con- 
struct LFs of each sample. However, we caution again that our 
Lir of AGN is total infrared emission of optically-selected AGN, 
and thus, it includes infrared emission from AGN host galaxies. In 
FigQ~3] we show LFs separately for AGN (red squares) and SFG 
(blue triangles). Composite galaxies are included in the AGN LF. 
The total IR LF shown in Fig[TT]is also shown for comparison in 
the gray dashed line. It is interesting that both AGN and SFG IR 
LFs have a break at around logLiR ~ 11.0. Then, the IR luminos- 
ity density of AGN exceeds that of SFG at logLiR > 11.0, with 
flatter slope toward the brigher end. In contrast, the faint-end slope 
is steeper for SFG. 

Next, we fit the double-power law (Eqs.[2]and[3} to the AGN 
and SF IR LFs, exactly as we did for the total IR LF in FigQj] The 
best-fit parameters are summarized in Table Q] As we mentioned 
above, the most notable difference is at the bright- and faint-end 
slopes, where the AGN and SFG LFs dominate, respectively. The 
L* of AGN is slightly brighter than that of SFG. 

One of the pioneers in IR AGN LF is 



iRush. Malkan and Spinoglid ( 1 19931) . who selected 893 IR 
galaxies (of which 118 are AGN) from IRAS 12/^m to minimize 
wavelength-dependent selection bias es towards bluer Seyfert 1 
nuclei and redder embedded AGN (Spinoglio & Ma lkanl 1 1989b . 
Despite small differences such as selection, AGN classification, 
and Lir based on the IRAS, their total IR LF of AGN in their 
Fig.9 is in very good agreement with ours, showing the same AGN 
dominance at larger Lir. 

What do these differences in LFs bring to the IR luminosity 
density, Qir, by AGN and SFG? We estimate the total infrared lu- 
minosity density by integrating the LFs weighted by the luminosity, 
separately for AGN and SFG. We used the double power law out- 
side the luminosity range in which we have data, to obtain estimates 
of the total infrared luminosity density, Qir, for AGN and SFG. 

The resulting total luminosity density (Qir) is, 



iL IR 



5.9+11 x 10 7 L Q Mpc- 3 , and Clfg N =1.9±°-* x 10 7 



1.3 ~ iW -£J B ±V±l^ , ilJJ -l.V_Q 7 

L Mpc~ 3 , as summarized in Table Q] Errors are estimated by 
varying the fit within 1 a of uncertainty in LFs. The results show 
that among the total IR luminosity density integrated over all 

the IR luminosity range, 75% (- 



i SFG 

aon'IoHpo ) of IR luminosity 

IR + ll IR 



density is emitted by the SFG, and only 24% ( „AGN R n SF6 ) is by 

n iR +n iR 

AGN at z=0.031. The AGN contribution is larger than z=0.0081, 
however, is still small compared with higher r edshift results 
jPerez-Gonzalez et alllgooll iLe Floc'h etaLll2005l ; iMagnelli et ail 
2009; Goto et al. 201(3). 



Once we have Qf^ G , we can estimate star formation den- 
sity emitted in infrared light. The SFR and Lir. is related by the 
following equati on for a Salpet er IMF, <f> (m) oc m -2,35 between 
0.1 - 100M Q dKennicuttlll998h . 



SFRlMeyr- 1 } = 1.72 x 10" 10 L/h[L g 



(4) 



By using this equation, we obtain SFR density = 1.3±0.2 

10" 2 MQt/r _1 . 

If we limit our integration to ULIRG luminosity range 
(Lir > 10 12 L Q ), we obtain, nf|' G (ULIRG)= 1.4±g;| x 10 5 
L Mpc -3 , and Qfg N (UURG)= 2.5tH x 10 5 L Mpc~ 3 . In 
other words, at ULIRG luminosity range, AGN explain 64% 

(if g N (ULIRG) 



(lfg N (U LI RG)+O.S^ G (ULIRG) 



) of IR luminosity, again showing 



the AGN dominance at the bright-end. 

Note that Qfg N (ULIRG) larger than Qjg al (ULIRG) is ap- 
parently inconsistent. This is due to the limitation in extrapolat- 
ing the double-power law to obtain Qir. For example, the AGN 
LF has a flatter bright-end slope of beta=3.1±0.2 than beta of 
3.54±0.09 for total LF. Once extrapolated to larger luminosity at 
logLiR > 12.5, the flatter slope could produce larger Qir once 
integrated even if all the observed data points for AGN are be- 
low those of total IR. To obtain more accurate estimate of Qir we 
clearly need data points in larger luminosity range than the AKARI 
observed. In this paper, we tried to obtain best-estimates of each 
Qir for AGN, SFG, and total. As a result, we caution readers that 
QJr" 1 is not necessarily the algebraical sum of Qf^ N and Q^r . 
In the LIRG luminosity (Lir > 10 Lq), results are 



Q S ir G (URG)= 2.1 



+0.7 



x 10° 



H0.5 



:)5 x 10 6 L Mpc 3 . This shows AGN contribution is already 



L Mpc -3 , and ^"(LIRG)= 



1.9 

down to 47% ( 
the LIRG range. 



if g N (LIRG) 



n ?H ™ (U LI RG)+Cl^ lJ (LI RG) 



) of IR luminosity in 
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Figure 14. Evolution of IR luminosity density by star-forming galaxies. Results from this work is plotted with stars at z=0.031. The red, blue and orange stars 
show IR luminosity density from all galaxies, from LIRG only, and from ULI RG only. Highe r redshift results in the solid lines are from the AKARI NEP deep 
field <Goto et al.l2010al) . Results at z=0.0082 are from the AKARI RBGS iGoto et all2010d) . Shown with different colors are IR luminosity density based on 
IR LFs (red circles), 8/jm LFs (stars), and I2fim LFs (filled triangles). The blue open s quares and orange filled squares are for only L I RG and ULIRGs, also 
based on our Ljr LFs. O verpl otted dot-dashed lines are estimates from the literature: Le Floc' h et al. (2005), Masnelli et al. (2009), Perez-Gonzale z et alj 
l2005t).ICaputi et alj j2007l). a nd lBabbedge et alj l200d) are in cyan, yellow, green, navy, and pink, re spectively. The purple dash -dotted line shows UV estimate 
bv lSchiminovich et alj l2005). The pink dashed line shows the total estimate of IR (IR LF) and UV i Schiminovich et al. 20051) . 



INFRARED LUMINOSITY DENSITY AND ITS 
EVOLUTION 



4.1 



Evolution of ilfj[ G 



We have separated the SIjr G from Qf^ N . Now we are ready to 
examine the evolution of Q,sfr without contribution from AGN. 
In Fig(T4l we plot the evolution of Qt r R as a function o f redshift. 
Higher redshift results are taken from lGoto et al who also 

tried to exclude AGN using SED fitting. Results from the Spitzer 
survey and GALEX survey are also plotted. The 0,f^ G shows a 
strong evolution as a function of redshift. The best-fit linear rela- 
tion is Qijj G oc(l+z) 41±0 ' 4 . We have removed Lir of optically- 
classified AGN in this analysis. A caveat in this process, however, 
is that Lm of optically-selected AGN perhaps includes some IR 
emission from star-formation in AGN host galaxies. In this sense, 
^ir G estimated here is a lower-limit. However, this should not af- 
fect nfn G more than 30% since the Q,f§ N is smaller than Qf£ G 
by a factor of 3 as we show in the next section. 

In comparison, our results are in good agreement with pre- 



vious works shown in th e dash-dotted l i nes in Fig[l4] In the 
form of Q. i R oc (1+z) 7 . iLe Floc'h et alj J2005l) obtained 7 = 
3.9 ± 0.4. Iperez-Gonzalez et al.l d2005l) obtained 7 = 4.0 ± 0.2 
iBabbedge et al. (2006?) obtained 7 = 4.5 + ° g. iMagnelli et all 



(2009) obtained 7 = 3.6 ± 0.4. 

Once the IR luminosity density is separated into ULIRG and 
LIRG contribution, we found Slfl a (ULIRG) oc(l+z) 10 ' 0±0 "\ 
and nf£ G (LIRG) oc(l+z) 6 - 5±0 ' 5 . fifJT (ULIRG) shows more 
rapid evolution than flf^ (LIRG), showing importance of lumi- 
nous IR sources at high redshift. 



4.2 Evolution of Qfg N 



In turn, we can also investigate fl'iR . By integrating IR LFagat 
inFig[l3] we obtained D.fg N =(l.9t° 5 T ) x 10 7 L Mpc~ 3 . 

In Fig[l5] we show the evolution of flj G , which shows 
a strong evolution with increasing redshift. At a first glance, 
both ilfji N and f2f^ G show rapid evolution, suggesting that 
the correlation between star formation and black hole accretion 
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Figure 15. Evolution of IR luminosity density by AGN. Results from this work is plotted with stars at the median redshift of z=0.031. The red, blue and 
orange point s show IR luminosity density from all AGN, from LIRG AGN only, and from ULIRG AGN only. Higher redshift results are from the AKARI NEP 
deep fi eld jGoto et al. 2010a), with contribution from sta r forming galaxies removed. Brown triangles are Q^GN com p U t e d from the 8fim LFs iGoto et al] 
l2010j) . Results at z=0.0082 are from the AKARI RBGS <Goto et alBoiOcb . 



rate continues to hold at higher redshifts, i.e., galaxies and black 
holes seem to be evolving hand in hand. When we fit the evo- 
lution with (1+z) 7 , we found ilfg N oc(l+z) 4 1±0 5 . A caveat, 
however, is that flf§ N estimated in this work likely to in- 
clude IR emission from host galaxies of AGN, although in opti- 
cal the AGN component dominates. Therefore, the final conclu- 
sion must be drawn from a multi-component fit based on better 
sampling in FIR by Herschel or SPICA, to separate AGN/SFG 
contribution to Lir. The contribution by ULIRGs quickly in- 
creases toward higher redshift; By z=1.5, it exceeds that from 
LIRGs. Indeed, we found Slfg N (ULIRG) oc(l+z) 8 7±0 ei and 
nfg N (LIRG) oc(l+z) 5 4±0 - 5 . 



5 SUMMARY 

We have cross-correlated the AKARI IR all sky survey with 
the SDSS to find 2357 IR galaxies with optical spectra. Using 
AKARI's 6-band IR photometry in 9, 18, 65, 90, 140, and 160/nm, 
we have measured Lir via SED model fitting. The AKARI's 6 
bands cover the crucial far-IR wavelengths across the peak of the 
dust emission, providing us with more accurate Lir measurements 



than IRAS. By separating SFG/AGN using optical line ratios, we 
constructed local IR LFs separately for SFG and AGN. We also 
computed local infrared luminosity density through the derived 
LFs, and compared 0.jr G and 0,fg N to those at higher redshifts. 
Our findings are as follows. 

• The local IR LF with the AKARI data at the median redshift 
of z=0.031 agrees well with that at z=0.0081, and shows smooth 
and continuous evolution toward higher redshift results from the 
AKARI NEP deep field. 

• By integrating the IR LF weighted by Lir, we obtain the local 
cosmic IR luminosity density of Qir= (3.8^2) x 10 8 L Mpc -3 . 

• The fraction of AGN and composite galaxies show continuous 
incrase from 25% to 90% at 9 < logLiR < 12.5. 

• The SFR based on self-absorption, and extinction corrected 
Ha correlates with Lir for SFGs. 

• H a [Hp ratio shows a weak increase with Lir, however with 
a large scatter. 

• The Ljoiii] correlates well with Lir for AGN. 

• The AGN contribution to Qir becomes dominant at Lir > 
1O U L , above the break of the both SFG and AGN IR LFs. At 
Lir sC 10 n L Q , SFG dominates IR LFs. 

• LIRG and ULIRG contribute to Q m a little; Only 1.1±0.1% 
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of Qir is produced by LIRG (Lir > lO 11 /^), and only 
0.03±0.01% is by ULIRG (L IR > 10 12 L Q ) in the local Universe. 

• Compared with high redshift results from the AKARI 
NEP deep survey, we observed a strong evolution of 

fi?£ G oc(l+z) 41±0 - 4 , and nfg N oc(l+z) 41±0 - 5 . 
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